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The spectral problem of the Dirac equation in an external quadratic vector potential is considered 
using the methods of the perturbation theory. The problem is singular and the perturbation series 
is asymptotic, so that the methods for dealing with divergent series must be used. Among these, the 
Distributional Borel Sum appears to be the most well suited tool to give answers and to describe 
the spectral properties of the system. A detailed investigation is made in one and in three space 
dimensions with a central potential. We present numerical results for the Dirac equation in one 
space dimension: these are obtained by determining the perturbation expansion and using the Pade 
approximants for calculating the distributional Borel transform. A complete agreement is found 
with previous non-perturbative results obtained by the numerical solution of the singular boundary 
value problem and the determination of the density of the states from the continuous spectrum. 
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INTRODUCTION. 



It has been argued since a long time [l| that the Dirac 
equation with an unbounded potential in vector cou- 
pling has no discrete but only completely continuous 
spectrum. It was shown in .2, 3] and subsequently con- 
firmed in literature, that precise conditions have to be 
imposed to the potential for this to be true. In two later 
papers, [3,[5[, Titchmarsh proved that the Dirac equa- 
tion with a linear vector potential satisfies the needed 
requirements and he studied the relativistic quantum 
mechanics of an electron in a constant (or piecewise con- 
stant) electric field. The first order of the asymptotic 
perturbation expansion was explicitly calculated due to 
the possibility of integrating the Dirac equation in one 
space dimension by means of hypergeometric functions, 
but hard technical difficulties prevented further analyt- 
ical developments, as well as the extension to a higher 
power law for the interaction or to more general inho- 
mogeneous electric fields. The Dirac equation with a 
vector potential that grows sufficiently fast at infinity 
is an incomplete dynamical problem. In classical terms, 
the particle can arrive at infinity in a finite time; in the 
operator language, we have (2,2) deficiency indices and 
the appropriate asymptotic boundary conditions must 
be defined. There is an infinite number of such bound- 
ary conditions of the self-adjoint type not very mean- 
ingful from a physical point of view. The most physical 
condition is the absence of sources at infinity, corre- 
sponding to the Gamow-Siegert conditions, that assume 
outgoing wave functions only: in this case, however, the 
Hamiltonian is not self-adjoint, the spectrum has com- 



plex eigenvalues, the dynamics is dissipative and the 
eigenfunctions have the meaning of metastable states 
of the dissipative dynamics itself: 
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where T„ = —2lm(E n /H) > 0, is the inverse of the 
mean lifetime of the metastable state ip n . Actually such 
states are stationary, if we neglect the decrease of the 
norm, so that they should be called metastationary; we 
prefer not to call them resonances, since they are indeed 
different from the usual resonances and have a direct 
physical meaning. 

Thinking of the transition from the non-relativistic 
to the relativistic system as a perturbation process in 
terms of the small parameter (1/c), the disappearance 
of the Schrodinger bound states connected with the con- 
fining potential defines a singular perturbation problem. 
A numerical non-perturbative investigation of the one 
dimensional Dirac equation with linear and quadratic 
potential has recently appeared Q : the purpose of that 
paper was to follow very closely the change of the spec- 
trum when passing to the relativistic regime and to de- 
scribe the spectral concentration @, [|| at finite values of 
(1/c), or, in a more physical language, to determine the 
density of the states of the relativistic system. It was 
indeed produced a numerical evidence that the spec- 
trum is completely continuous and that it is given by a 
sum of Breit-Wigner lines reducing to 5-functions cen- 
tered at the non-relativistic eigenvalues for smaller and 
smaller values of the ratio of the interaction to the rest 
energy, so that the spectral measure becomes atomic 
as it should. In Q the pair production rate was also 
calculated from the line width of the lowest state find- 
ing a perfect agreement, for the linear potential and 
in the range of validity of the first perturbation order, 
with the results obtained from the imaginary part of 
the Schwinger effective action, [9[ . It should also be no- 
ticed that the QED results on the pair production for 



fermions in non constant electric fields are not yet so 
sound, although interesting proposals can be found in 
literature, (see [10( for an up to date review). We can 
finally remark that relativistic models with power-law 
potentials, in three space dimensions and in spherical 
geometry, have been used in the study of composite sys- 
tems like Quarkonium, in order to determine the mass 
spectrum of some meson families, [llj ]. 

In this paper we come back to the spectral properties 
of the Dirac equation with a quadratic vector poten- 
tial in a perturbation framework, completely different 
from that adopted in Q . Indeed, since the quantitative 
and numerical results on this subject are rather new, 
we find it interesting to make a comparison of indepen- 
dent computational approaches and eventually to have 
a confirmation of the results. In particular, we use here 
the method of the Distributional Borel Sum (hereafter 
DBS) , that has proved a very useful tool for dealing with 
physical systems of singular nature. Originating from a 
suggestion given by 't Hooft for double well problems, 
[ig(, the DBS was studied in a series of papers [T^|-[TH 
and it was successfully applied to A$ 4 lattice fields with 
large coupling in (l2| , to the non-relativistic Stark res- 
onances in |1 4| an d to the double well Schrodinger op- 
erators in [131 Il5|. To our knowledge however, despite 
the very deep mathematical development of the sub- 
ject, no explicit numerical calculations with the DBS 
have been made up to date and no applications of it to 
the Dirac equation have been considered. In this sense 
this paper represents a novelty and a test of the method 
for its possible practical uses. We will also show that 
the basic results of [l4[ can be brought to bear to our 
present context, although some peculiarities due to the 
relativistic nature of the problem will emerge and must 
be taken into account. 

The plan of the paper is the following. In the next 
section, for the sake of completeness and because the 
knowledge of the DBS is not so widely diffused, we 
give a brief description of the method and we state 
more properly the spectral problem. In section IIIII we 
study the general conditions for applying the DBS to 
the Dirac equation and we discuss the possible strate- 
gies for the one-dimensional and the three-dimensional 
problems. Finally, in section HVl we present the numer- 
ical treatment for the one-dimensional case. Using a 
perturbation parameter proportional to the ratio of the 
interaction to the rest mass energy, we calculate the 
perturbation series of the fundamental state energy up 
to a large order and we then construct its Borel trans- 
form determining its asymptotic behavior. We then ap- 
proximate the Borel transform by Pade approximants. 
Although, in a strict mathematical sense, this is not 
a completely rigorous procedure since it is not proved 
that the Borel transform is a Stieltjes function, how- 
ever an idea of the accuracy of the approximation can 
be obtained by comparing the position of the poles of 
the Pade approximants with the location suggested by 
the asymptotic form of the perturbation series: indeed, 
as shown in section [TV] we find the Pade poles exactly 
where the asymptotic behavior indicates they should be. 



We finally check the stability of the poles with the in- 
creasing orders of the Pade approximants, implying the 
stability of the imaginary part of the perturbed energy. 
The latter is then calculated and it is found a perfect 
agreement with the results of Q, confirming them and 
proving the correctness and the effectiveness of both 
methods also in explicit numerical calculations: this is 
even more interesting in view of the fact that each of 
the two methods is more efficient in a different range 
of the values of the coupling constant, so to allow for 
the choice of the most appropriate one in a practical 
situation. 



II. THE SPECTRAL PROBLEM. 



We consider the (3+l)-dim Dirac equation interact- 
ing by means of a central vector potential (0, V(r)). It 
is well known that the use of the spherical spinors, [lj|, 
easily leads to the diagonalization of the angular mo- 
mentum and to the reduction of the Dirac equation to 
the system of the two first order differential equations 

-(W + mc 2 - V(rf)Xi(r) - h(jj- + -)x 2 (r) = 
h ( 4- - ~)Xi(r) + -(w-mc 2 - V(r))x a (r) = 

(III) 

In (jll.l[) to is the mass and W the energy of the par- 
ticle, while X\(r) = rg(r) and -^(r) = rf(r), where 
f(r), g(r) are respectively the 'large' and the 'small' 
spherical component of the spinor. Finally the param- 
eter k accounts for the angular momentum and the 
parity: k = — it + 1) for j = £ + 1/2 and n = £ for 
j = £ — 1/2, |16j, so that, in three space dimensions, k 
can assume all integer values except zero. If, instead, 
we let n — in equation (jll.l[) and we change r in x, with 
— oo < x < +oo, we obtain the one space-dimensional 
Dirac equation that has been discussed in @. We now 
assume a quadratic potential 

V{r) = (1/2) mwV, 

we rescale the system by introducing the dimensionless 
variables 

/"muj\ 1 / 2 1 / frw x 1 / 2 

V = T \— ) 0= 2^) 

E = ^-{W -mc 2 ) 
nu> 

and we define the unknown functions 

$ 1 (r)=2- 1 /2( Xl (r)+7X 2 (r)), 
$ 2 (r) = -i2- 1 / 2 (X 1 (r)-iX 2 (r)) 
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The system pi.l[) becomes then 
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dy 



<I> l(/y ,-if 7 l+ / 3(S-j / 2 )) 
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^ + /3(^-y 2 ))<f 2 (y) = 

By eliminating $2 (y) we find the second order equation 
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K l + 2(3 2 {E~y 2 ) 
y {y-2inP) 
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$1(2/) =0 



(II.2) 



Besides infinity, this equation, with y S [0, +00), 
presents the obvious singularity at the origin, absent in 
the one-dimensional case with k = and y £ (—00, +00). 
The non-relativistic limit for /5 — >■ 0, 
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E-y- 
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=0 (II.3) 



reproduces the usual Schrodinger equation, with orbital 
angular momentum £, in a quadratic potential. The 
properties of the equation with [3 7^ 0, however, are 
very different from those of the non-relativistic coun- 
terpart, since pi.2[) . even neglecting the presence of an 
imaginary term, presents an asymptotic oscillatory be- 
havior that prevents the existence of any normalizable 
solution. Therefore, as we said, the perturbation ex- 
pansion from the non-relativistic system turns out to 
be singular. 

We can transform the domain of definition of the dif- 
ferential equation (|II.2[) with a map U defined by 



$(tf) = l/(*i(»)) 



1 



*i(y)- (n.4) 



E _ y2 _ K (n + l) 



The transformed differential equation becomes then 

d 2 
dy 

V (y-2iK(3)y i 
2n{y 2 -E) 3 k 2 



k(k + 1) 



V 



;) +f 3 2 (( y 2 ~E) 2 + 



(y-2in(3)y (y - 2 i k (3) 2 y 2 



Hy) = 

(II.5) 



Observing that (3 > 0, we find it useful to introduce 
the parameter g = =fi{3 according to whether k > or 
k < respectively. With such a choice of the signs, 
when g > no additional singularities besides the origin 
and the infinity appear at finite values of y and the 



transformation U is unitary. The equation to be studied 
can finally be written as 



H g $(y) = 0, H g = --^ + V(E,g,K,y) (11.6) 



whose 'potential' is given by 

V(E, g, k, y) = V (E, K , y) + Vi (E, k, g, y) (IL7) 
where 

2 H(M±i) 

Vo = y --EH 72 1 



Vi =±2 5 1/ 
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N(N±i) 
(y + 2M<?)y 2 



/ .,.2 ^2 ± ^\n\(y 2 -E) 3|«| 2 

(y + 2|At|fl)i/ (y + 2| K | 5 ) 2 ?/ 2 
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The double signs correspond again to k > and k < 
respectively, while the dimensionless energy E, entering 
V(E, g, k, y) in a polynomial expression, can be consid- 
ered, for the moment, a parameter that will be deter- 
mined during the calculation, ft is common sense - 
and it can also be proved rigorously - that the equa- 
tion (|II.6|) with positive g admits bound states: start- 
ing from these we can then try to make an analytic 
continuation in the complex plane of the coupling con- 
stant g from the real positive to the positive or nega- 
tive imaginary axis in order to recover the values of the 
parameters of the initial problem. The procedure, how- 
ever, presents some delicate points. Indeed the analytic 
continuation of the eigenvalue of the complete equation 
would be immediate if the eigenvalue itself was given as 
a convergent series expansion in the coupling constant 
g. Unfortunately this is not the case: the perturbation 
expansion in g is asymptotic and has zero radius of con- 
Since finally our original problem has 



21] 



vergence, 

exactly y 2 <0, the system we are dealing with is there- 
fore close to an unstable quartic oscillator, but for some 
non polynomial terms for which an additional discussion 
is in order: all the considerations we developed in the 
Introduction about the incompleteness, self-adjoint ex- 
tensions and boundary conditions at infinity thus apply. 

Let us now recall that the study of the asymptotic ex- 
pansions for treating perturbation problems in quantum 
mechanics was systematically undertaken since the be- 
ginning of the seventies and the main concept allowing 
to deal with divergent series was the Borel summability 
@, [LlL Ull . Given a formal series 



n=Q 

we define its Borel transform as 

00 



(II. 



(II.9) 



n=0 



where v is a fixed parameter independent of n that is 
usually chosen in relation to the asymptotic behavior of 
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the coefficients {o„}, but whose choice is otherwise irrel- 
evant. If then B(u) has a nonvanishing radius of con- 
vergence, admits an analytic continuation to a neigh- 
borhood of the positive real axis and if the integral 

/•oo 

f{z) = z- 1 / du {u/zf exp(-u/z) B(u) (11.10) 
Jo 

converges in a so called "Nevanlinna disk" Cr = 
{Re (z^ 1 ) > J2 _1 , R > 0}, we say that f(z) is the sum 
of S in Cr. In some favorable cases, where stronger 
properties hold, more convenient summation methods 
have also been proposed as, for example, the Stielt- 
jes sum, 18]: this can be defined when the classical 
Stieltjes moment problem has a solution and the result 
can be conveniently expressed in terms of a convergent 
Pade approximation. There are many cases, however, in 
which the Borel summability and, a fortiori, the Stielt- 
jes summability cannot be applied: notably, this occurs 
when the Borel transform develops some singular point 
along the positive real axis, so that the integral (|II.10[) is 
not defined. The DBS originates from the need of find- 
ing a more flexible method of summation that could be 
able to avoid the difficulties connected with the emer- 
gence of the singularities. The idea is to follow closely 
the way used in the discussion of the Riemann-Hilbert 
boundary value problem with data assigned on the pos- 
itive real axis. More precisely, we reduce the require- 
ment of the analytic continuation of B{u) in a neigh- 
borhood of the positive real axis to the assumption of 
the existence of the analytic continuation of B(u) in the 
intersection of a neighborhood of the positive real axis 
with the upper complex half-plane, so that the bound- 
ary value B(u + iO) is well defined for any u > 0. The 
integral yielding f(z) is then replaced by 

POO 

(P( z ) = z- 1 / B(u + i0) p(u/z)du, (11.11) 
Jo 

where the measure p(u)du is finite and has positive 
moments pk = / °° u k p(u) du. The function (f>(z) de- 
fined by the integral (jll.lip is analytic in the upper 
half-plane and it is called the upper sum. In particular 
for fj,k = r(fc + l + we have p(u) = u"exp(— u) and 
(III. lip reproduces the relationship between boundary 
values of analytic functions and many distributions de- 
fined by the inverse Laplace transform, (20| . Since then 
B(u — iO) = B(u + iO) we also define the lower sum 

/•OO 

4>(z) — z^ 1 / B(u + iO) p{u/z) du . 
Jo 

and we consider the real and the imaginary parts of 

<t>{z): 

d{z) = ^(<l>{z)-W))=^{<t>{z))- (H.12) 

It is then natural to assume f(z) as the DBS of the se- 
ries, while the discontinuity along the positive real axis, 



d(z), is uniquely determined and it has a zero asymp- 
totic power series expansion. In [23j d{z) was related to 
the Schwinger effective action. 



III. APPLICATION OF THE DBS. 

In this section we study the large order perturbation 
approach to the equation (|II.6|) in one and in three space 
dimensions and we stress some physical relations be- 
tween angular momentum and perturbation expansion. 
For the sake of completeness we will recall some useful 
facts related to general and well known properties of 
operators: greater details can be found in 0, H, HE HH • 
These notions can be collected in two different main 
subjects. The first concerns the analytic properties of 
the operator family, related to the coupling constant of 
the perturbation part. This is essential if we want to 
explore the changes of the spectrum of the operators 
during the analytic continuation process. The second 
one deals with the estimates of the asymptotic growth 
of the perturbation series, that imply the existence of 
appropriate summation mechanisms. 

The differential equation defined in pi.6|) and depen- 
dent upon a parameter g, has nice properties for g = 0, 
in the sense that we are able to determine exactly its 
discrete spectrum. The natural question to be posed is 
whether an eigenvalue p of Hq = H g= o with a certain 
(algebraic) multiplicity M gives rise to nearby eigen- 
values p\{g), p s (g) with total multiplicity M when 
the perturbation is switched on. If this is the case, 
the eigenvalue p is said to be stable. We just recall 
that the discussion of this type of problems is simple 
enough when the perturbation potential V\(E, re, g,y) 
is ifo-bounded, namely when the domain of defini- 
tion T>{V\) contains V(H ), since this property yields 
\\Vi4>\\ < a ll^o0|| + b\\4>\\t for constant a, b and for any 
4> € T>(Hq). From the previous inequality, as in the 
proof of the Kato-Rellich theorem, we can easily de- 
duce that T>(H g ) is independent of g and that, for any 
4> G T>(H g ), H g (j) is an entire analytic function of g. In 
this case the perturbation scheme is said to be regular. 
A slightly generalized strategy must be adopted when 
the perturbation is not i/o-bounded and specially when 
the behavior of the whole system depends dramatically 
upon the sign of the coupling constant or its square, 
as in our specific problem. The idea is to extract the 
fundamental properties of the relatively bounded case 
by means of the notion of analytic family of operators, 
introduced by Kato Q: this term denotes a collection 
of operators {Hg\ ge R dependent upon a coupling con- 
stant g taking values in a region R of the complex plane, 
where each H g has a non-empty resolvent and is de- 
fined on a domain T> independent of g; moreover, for 
each (j) G T>, H g 4> is required to be strongly analytic 
in g. Although weaker than relative boundedness, this 
analyticity property is still sufficient to study the con- 
tinuation of the eigenvalues in a region of the complex 
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plane of g and it can be proved that for £ in an ap- 
propriate open region of the resolvent set of H 9g the 
resolvent (£ — H g )~ 4 is analytic in g when \g — g \ is 
sufficiently small. The stability of discrete and non- 
degenerate eigenvalues immediately follows. For iso- 
lated eigenvalues of multiplicity n there exist t families 
of eigenvalues \i(g) admitting Puiseux expansions in 
terms of g 4 ^ qi , for integers {<ft}i=i,i) with total multi- 
plicity n. We are therefore left with the need to prove 
that a family of operators is analytic: the technique 
commonly used for this purpose is to establish appro- 
priate quadratic estimates that allow us to deduce that 
each operator of the family is closed and that the per- 
turbation term is small - in the sense of Kato 0] _ with 
respect to any member of the family, so that the do- 
mains are indeed independent of the coupling constant. 
This argument can be formalized as follows [8|: if To 
and Ti are closed operators with V{Tq) fl T>(T\) dense, 
then Tq + T± is closed if the inequality 



a = (1 - 2Eg 2 



-1/4 



(III.3) 



T + TiJ [T + TiJ > a[T* T + T* Tij + b 

(IIL1) 

is satisfied for some constants a, b. The one-dimensional 
quartic oscillator with a Hamiltonian H =p 2 +x 2 +g 2 x 
describes a situation very close to the Dirac equation 
in one space dimension we want to discuss and its spe- 
cific estimate, (25[, can be formulated as follows: for 
a < l-\Re{g 2 )\/\g 2 \ and for all 4> £ V{p 2 ) n V{x A ) 
there exists b such that a (\\{p 2 +x 2 )<j)\\ 2 + \g\ 4 \\x 4 (J)\\ 2 ) < 
\\(p 2 +x 2 +g 2 x 4 )<j)\\ 2 +b\\(l)\\ 2 . We finally observe that 
this estimate implies the Herglotz property of \i(g 2 ), 
namely Im(Ai(<? 2 )) > for Im(g 2 ) > 0, that gives pre- 
cise informations about the analytic structure of \{g 2 ) 
and shows that g = cannot be an isolated singular 
point, but it must be a limit point of singularities of 
the eigenvalues. 

We will now discuss the application of these general 
facts to the Dirac equation in a quadratic vector poten- 
tial. 

(a) The case in one space- dimension. 

Let us consider first the one-dimensional problem. If 
we put k = in pi.6l - III.7|) . we have the differential equa- 
tion 



The inverse relations are 



E = 2X (1 + Xa 2 ) 1 / 4 [ 1 + (1 + Xa 2 ) 1 ' 2 ] _1 



g= U -{l + Xa 2 fl\ 



(III.4) 



and by means of (|III.4|) it is possible to reconstruct the 
function E(g) from A(er). In the new parametrization 
P11.3p . equation (|III.2|) becomes 



2 u 4 

V + ay + —y 



$(y)=A$(y) (III.5) 



and the quadratic estimate has to be done for the opera- 
tor on the left hand side of (|TTL5|) . For <E V(p 2 )nV(y 4 ) 
and Im(er) ^ we consider the closed operators 



T Q =P 2 



Re(a)y , 



Ti = i1m(<r)y + —y 4 



(III.6) 



Everything is known about Tq since it is just the Hamil- 
tonian of a harmonic oscillator in the displaced variable 
(y + Re(cr)/2). A direct calculation given in Appendix 
A proves the quadratic estimate (jlll.ljt . 

It will be shown in Section [TVl that for parity reasons 
only even powers of the coupling constant a contribute 
to the perturbation expansion of the eigenvalues. As 
a consequence of (|A.1|) we can state that on the do- 
main V{p 2 ) n V(x 4 ) and with T , T x given in (|III.6|) . 
the operators H a = Tq + T\ form a holomorphic fam- 
ily with compact resolvents for a 2 in the complex plane 
cut along the negative axis. We thus get a perturbations 
series for each eigenvalue 



X(z) = 2J a n z r 



(III.7) 



n=0 



with z = a 2 . For (|III.7I) there exists a Borel transform 
(III . 9|) . The corresponding sum is given by 



X{z) = z- 1 -" J u v cxp(-u/z) B{u + iO) du (III. 



p 2 + (1 - 2Eg 2 )y 2 + 2gy + g 2 y 4 J $(y) 

= E(l - g 2 E) $(y) (III.2) 



defined for — oo < y < oo, where p = —i(d/dy). As 
we said, in order to apply the general theory, we need a 
quadratic estimate analogous to that of the anharmonic 
oscillator. A less cumbersome notation is obtained if we 
rescale the variables by a dilation y — > ay, p — > p/a and 
then we choose 



X = Ea 2 (l- Eg 2 ), 



2ga 3 , 



(b) The case in three space- dimensions. 

We next consider the three-dimensional problem as- 
suming k = — 1 in order to simplify the notation: indeed 
the treatment is qualitatively the same for any other 
allowed value of n. The equation (jII.6HII.7l) . then, spec- 
ifies to 



H a Hy) = 



2 +y 2 -2gy- 
3 

y 2 (y+2.9) 2 



9 2 {{y 2 -E) 



2(y 2 -E) 
y{y+tg) 

E$(y) (III.9) 
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and since g is proportional to c _1 the limit g — > cor- 
rectly reproduces the Schrodinger equation pi.3[) with 
angular momentum I = 0. 

With respect to (|III.2[) , equation (|III.9[) contains non 
polynomial terms proportional to g 2 . The rational 
terms have vanishing denominators when both g and 
y tend to zero and the simultaneous limit is not well 
defined. The expansion in g of the perturbation poten- 
tial leads to terms with higher and higher divergences 
in the origin: 



V x ~-2gy + g 2 [y- 4 (3 + 2 E y 2 + (E 2 - 2) y 4 
-2Ey e + y s ) - gy- 5 (i (3 + E y 2 - y 4 )) + 0{g 2 ) 



These terms make rather problematic the possibility of 
achieving a reasonable quadratic estimate. On the other 
hand the expansion in y in the neighborhood of the 
origin gives, up to to terms of 0{y), 



3 3 -4 Eg 2 9- 8 Eg 2 + 16 E 2 g 4 
V 



4 gy 



16 g 2 



(111.10) 



As the quadratic divergence at the origin is due to the 
angular momentum of the system, the term (3/4) y -2 
shows that, as a consequence of relativity, we are dealing 
with a spin (1/2) fermion and that the spin is the only 
angular momentum left when n = — 1. Besides these 
physical observations, however, the expansion cannot 
be used perturbativcly, not even for searching the local 
solutions at the origin. A different approach has there- 
fore to be searched. We thus find it necessary to fix 
the final value of the parameter g and to put this final 
value in the singular terms: we say that we freeze 
the parameter at its final value wherever it is needed. 
We then write the equation 



p + y + 



3 

4y 2 

-(y 2 ~E) 2 + 2 



y 



2gy + g 
2 E 



(y + 2g*)y 



3 y + 4g 
4<? 2 (y + 2g*) 2 y 

(111.11) 



that coincides with (|III.9[) when g* = g. We denote 
respectively by Tq and T± the operators in the first and 
in the second square bracket; we then observe that with 
g frozen at any non- vanishing g* the factors (y + 2g*)~ 1 
are bounded in < y < +oo and do not need further 
estimates. Therefore the two terms of T\ diverging like 
y~ 4 at the origin and presenting an at most constant 
behavior at infinity are small in the sense of Kato with 
respect to y 2 + (3/4) y~ 2 in T , so that they can be 
neglected in the estimate. Since it is well known that 
p 2 -\-y 2 + J (j + l)y~ 2 is a closed operator on the maximal 
domain of the functions satisfying the condition $(y) ~ 
at y = 0, that correspond to the regular solutions 
of the differential equation (|III.9|) - and therefore To 
is closed on the maximal domain with <&(?/) ~ y 3 ^ 2 at 
the origin - , by means of a calculation analogous to 



(jA.ip we conclude that To + T\ is closed and that the 
set {H g }, whose elements are specified in (IIII.9I) . form 
an analytic family of operators with compact resolvent. 
We can thus conclude that for the three-dimensional 
case also we get a perturbation series like (|III.7|) with 
z = a 2 for which the Borel transform (|II.9|) and its 
inverse (jIII.8l) are well defined. 

It is evident that the concrete calculation of the per- 
turbation expansion is much more difficult in three di- 
mensions than in one. In practice the expansion (|III.10p 
could suggest as preferable the use of a variational 
method by defining the isospectral dilated operator 



H n (g,E)= v 

V- 



p + y + 



l)y -<!) 

+.9V (v 2 y 2 - e) 2 



3 
4y 2 
3(3ryy + 4<7) 



^(j)y + 2g) 2 y 



Igrfy 



1 2 

TV 



(vy + 25)?/ 



where g = —i\g\, r\ = exp(iz9), < < 7r/6. We could 
then approximate the first eigenvalues restricting the 
Hamiltonian on the linear space of the first n eigenvec- 



tors of p 



y 



(3/4) y . We leave this problem for 



future investigations. 



IV. NUMERICAL DEVELOPMENTS. 

In this last section we present the numerical calcula- 
tions concerning the perturbation treatment of equation 
(|III.5|) yielding the determination of A(er) from which we 
can obtain E{g) by inverting the definitions (|III.3[) . We 
take p 2 + y 2 as unperturbed operator and we write the 
perturbation part in the form aUi+a 2 U2, where U\ = y 
and U2 = y 4 /4. We introduce the usual creation and 
destruction operators 

a = 2-Va (d/dy + y), a+ = 2" 1 / 2 (-d/dy + y) , 
[a, a 1 ] = 1 , 

so that 



p z + y 2 = 2a T a + 1, 



2-V3 (at+o), 



y4 = 4-I ( fl t 4 + 4a t 3 a + 6a t 2 a 2 + 4a t a 3 

+ a 4 + 6at 2 + 12at a + 6a 2 + 3). 

In the standard Dirac notation, the solutions of the 
eigenvalue equation for the unperturbed operator, 
(2cva + l)\ip) = A\ip), are the usual occupation num- 
ber states |n), expressed in terms of Hermite functions, 
where A = A„ = 2n+ 1 with integer n. For later use we 
recall the matrix elements of y, y 2 and y 4 with respect 
to pairs of such states, namely 

(k\ y \n) = 2~ 1/2 (Vn+l5 k . n+1 + y/nSk, n -i) 
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(% 4 |n) = 2- 2 (V(ri + 4)(n+3)(n + 2)(n+l)4,„+4+ 
(4n+6)V(»+2)(n+l)<5»,, n+ a + 
(6n 2 +3(2n+l)) 4,„ + (4n-2)Vn(n-l)4,n-2 + 
A /n(ra-l)(n-2)(n-3)4,n-4) (IV.l) 

where <5j,j is the usual Kronecker symbol, taking values 
one or zero according to whether the indices are equal or 
different. Since the Hermitc functions have the parity 
of n, we see from (|IV.1[) that (k\Ui\n) is vanishing when 
the integers k and n have the same parity and (k\U~2\n) 
is vanishing when the parity is opposite. We now set 
up the standard perturbation framework by defining 



A„ = (2n + l)+ Efl'A. 



(i) 



\lpv 



As usual, the global phase of \ip n ) can be chosen in such 

a way to have (n\ipn } = 0. A straightforward calcula- 
tion leads then to the first and second order quantities 

A« = 0, A<*> =2" 4 (6n(n+l)-l), 

|V£ X) ) = -(2V2)- 1 (V^+T |n+l) - |n-l)) , 



128 

, (2n+l)V(n+2)(n-l) | 



32 

(2n+l)^?i(n-l) 

32 

A/n(n-l)(n-2)(n-3) 
128 



-2) 



|n - 4) (IV.2) 



and to the recurrence relation, for i > 2, 

"(p 2 + y 2 )-(2n+l)] [7-^-1))+ 
i-l 

(IV.3) 

We now prove by induction that for any m > we 
have: 



(a) A 

(b) |V 



(2m+l) _ 



0. 



(2m+l)\ 



has the opposite parity of n, 



(c) |^n 2m+2 ' ) ) has the same parity of n. 



The relations PV.2P give the initial step of the inductive 
argument. Suppose now the properties we require are 
true for i < 2m. The scalar product of ()IV.3[1 by |n) 
gives 

A» = (nlW^) + (n\U2\^t 2) ) (IV.4) 

and we first consider i = 2m + 1. By the inductive 
hypothesis both the matrix elements of U± and {/ 2 are 



(s) 

vanishing. The last sum is also vanishing, since A« = 
for odd s and {n\ipn m+1 s ') = for even s. As a conse- 
quence Xn™ 1 ^ — 0. Looking then at the perturbation 
contributions to the states, we have 



IV#>> = I>>W-- 



(IV.5) 



fc=0 



where the scalar product of PV.3P by \k) provides the 
relation, for k =/= n, 

(fc|4 l) ) - (2n-2fc)- 1 ((fc|C/ 1 |^ ( r 1) ) + 

(k\U 2 \^t 2) ) - E Ai s) (£#<T s) )) (IV.6) 

8=1 ' 

from which the parity properties are easily deduced. 
The induction is therefore complete and it implies, in 
particular, that the perturbation expansion of the eigen- 
values is in <t 2 rather that in a. 

By using (jIV. HIIVT6]) we can set up a recursion scheme 
by which the expansions of the eigenvalues are deter- 
mined. With the help of some computer algebra we 
have calculated the first 94 coefficients of the pertur- 
bation expansion in a of the lowest unperturbed eigen- 
value and of the fundamental state in terms of rational 
numbers, namely with infinite precision. Starting from 
there we have continued the expansion with a fortran 
code in floating point up to 250th order, thus obtaining 
a contribution to the eigenvalue up to order 125 in a 2 . 
The coefficients of the asymptotic series in a 2 have al- 
ternate signs and, on the basis of the numbers we have 
calculated, we find the evidence for an asymptotic be- 
havior 



a„«|(l-^)(-ir(|)"r(r, + l/2) (IV.7) 

To be more precise, we have defined the sequence {a n } 
where a n is given by the ratio of the coefficient a n of 
the initial divergent series divided by the right hand 
side of equation pV.7|) . We have then used the Shanks 
transformation A n = (a n+ ia n -i - a 2 n )/(a 

n+l "T (In— 1 

2o n ), [271] . in order to improve the convergence of {a n } 
and we have found A125 = 1.00009, where the rounding 
numerical error can be estimated of the order of 10 -4 . 

Some comments on the general features of the per- 
turbation series and of its Borel transform are in order. 
For positive values of a 2 the series is Borel summable. 
For negative a 2 (namely for imaginary a, which is the 
case we are interested in), all the terms of the series 
acquire equal signs and the series itself becomes Borel 
summable only in the distributional sense. We will 
therefore assume the asymptotic series with all positive 
terms and thus consider z = a 2 = 4 f2* > in the Borel 
anti-transformation integral (|III.8[) : the meaning of £7» 
will be clarified below. Dropping (—1)" in the asymp- 
totic behavior pV.7[) it is clear that the Borel transform 
develops singularities on the positive real axis, the first 
of which has to be expected in the neighborhood of 8/3. 
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The shift in the energy of the state, given by f(z) in 
(III.12[) . can be plainly calculated by taking the princi- 
pal part of (|III.8P with respect to the poles on the real 
axis. 

The imaginary part of the perturbed energy, given 
by d(z) in (|II.12[) , has a much more interesting physi- 
cal meaning: indeed 2d(z) gives the decay constant T 
of the state itself as in and, in an interpretation 
connected with a second quantization framework, to 
the particle-antiparticle production rate iv? — 2Im£ c ff, 
where £ e ff is the effective Lagrangian of the system, 
[H, Q . The imaginary part will be computed by sum- 
ming the contributions of the positive real poles to the 
integral pil.8|) , calculated along the integration path 
encircling those poles in the upper complex half plane. 
An efficient way to make the calculation is to use a Pade 
approximation for the Borel transform B(u) of the per- 
turbation series. Since B{u) is expressed by a series 
expansion up to order 125, we can take Pade approxi- 
mants of rather high order and to look at the stabiliza- 
tion of the values of the poles. The results for the seven 
lowest poles are summarized here below. 



[58,59] 


[59,58] 


[59,59] 


[59,60] 


[60,59] 


[60,60] 


2 


67561 


2 


67562 


2 


67571 


2 


67575 


2 


67576 


2 


67572 


2 


71528 


2 


71529 


2 


71583 


2 


71611 


2 


71612 


2 


71589 


2 


79311 


2 


79314 


2 


79480 


2 


79563 


2 


79564 


2 


79497 


2 


92226 


2 


92235 


2 


92678 


2 


92899 


2 


92901 


2 


92722 


3 


12725 


3 


12752 


3 


13965 


3 


14582 


3 


14589 


3 


14082 


3 


45803 


3 


45883 


3 


49492 


3 


51474 


3 


51500 


3 


49830 


4 


03675 


4 


03958 


4 


17104 


4 


26976 


4 


27146 


4 


18431 





[57,57] 


[58,58] 


[59,59] 


[60,60] 





.05 


.251588e- 


5 


.251584e- 


5 


.251588e- 


5 


.251586e- 


5 





.08 


.321991e- 


3 


.321360e- 


3 


.321933e- 


3 


.322124e- 


3 





.10 


.162741e- 


2 


.161602e- 


2 


.162615e- 


2 


.163092e- 


2 





.12 


.483648e- 


1 


.476170e- 


1 


.482741e- 


1 


.486376e- 


1 





.15 


.146167e- 


1 


.141402e- 


1 


. 145544e- 


1 


.148164e- 


1 





.20 


.456068e- 


1 


.426354e- 


1 


.451921e- 


1 


.470098e- 


1 





.25 


.925871e- 


1 


.837260e- 


1 


.913049e- 


1 


.970549e- 


1 



where, for different values of f2, we have reported the 
numerical data for the imaginary part of (|III.8|) , namely 



Im(A) 



(40 a ,) 1 + y 



E 



Res P AP 



{poles Pi} 



(IV.8) 



with v = —1/2. In PV.8P B(u) is the Borel transform 
(III.9[) of the perturbation series and P[ m ,n] is the [m, n] 
Pade approximant of B(u) with poles {pi}. By Res Pi 
we have indicated, as usual, the residue at the pole pt. 
Finally the value of ft* has been constructed by using 
the formulas (|IlL3|) with ft* = ft (1 + 2Efl)' 3 / 2 and 
ft = \g\ 2 is the same parameter used in Q. We strongly 
stress that while in (|III.3|) . with imaginary g and a, 
A = E to the order |<7| 2 , the relationship between a and 
g is not so straightforward, as the correction due to the 
term a 6 is large and cannot by no means be neglected: 
on the contrary, it proves fundamental to establish the 
correct relationship between the perturbation data and 
the data calculated in Q , taking the zeroth order A = 1 
inside a. 

We also report a further check of stability of the 
results: the numerical values give the imaginary part 
pV.8p calculated with the use of the [60, 60] Pade ap- 
proximant for different values of the constant v. 



From these numbers it appears, first of all, that the 
locations of the lowest poles are in the neighborhood of 
8/3, as deduced from pV.7|) : this was to be expected 
and indicates the consistency of the approximation, as 
we said in the Introduction. We then see that their val- 
ues tend to stabilize with the increasing order of the 
Pade, the lower poles obviously faster than the upper 
ones, for which approximants of higher order would be 
required. Therefore, recalling that no further informa- 
tion can be deduced by the Pade construction when 
the sum of the degrees of numerator and denominator 
exceeds the order of the series, a considerably larger 
number of perturbation terms should be calculated in 
order to stabilize the higher poles. This means, in turn, 
that the calculations of the real and imaginary parts 
of the perturbed energy give much more precise results 
for small values of z, for which the contribution of the 
lowest poles is largely dominant: a natural fact in a 
perturbation framework, that can be clearly observed 
in the next table, 





i/=-l/2 


v=Q 


^=1/2 


v=l 





.05 


.251586e- 


5 


.251586e- 


5 


.251586e- 


5 


.251586e- 


5 





.08 


.322124e- 


3 


.322124e- 


3 


.322099e- 


3 


.322113e- 


3 





.10 


.163092e- 


2 


.163309e- 


2 


.163333e- 


2 


.163335e- 


2 





.12 


.486376e- 


1 


.488771e- 


1 


.489565e- 


1 


.489141e- 


1 





.15 


.148164e- 


1 


.150256e- 


1 


.151061e- 


1 


.150456e- 


1 





.20 


.470098e- 


1 


.485809e- 


1 


.487832e- 


1 


.484381e- 


1 





.25 


.970549e- 


1 


.101983 




. 100359 




. 100484 





We will conclude the section by comparing the results 
here obtained with those presented in [6j . The methods 
used to get them, as we said, are completely different. 
On the one hand, we have a perturbation approach that 
becomes more and more precise when the perturbation 
parameter decreases and needs many further terms for 
intermediate values of that parameter: the imaginary 
part of the energy is directly calculated by using the 
Pade approximants to invert the Borel transform in the 
distributional sense and integrating along the appropri- 
ate path in the complex plane. On the other, the deter- 
mination of the spectral concentration is obtained by 
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0,1 
Im(Ji) 
0,05 



O Data from Breit-Wigner width, [6] 
Data from Perturbation Theory 




0,15 0,2 



0,25 



FIG. 1: Imaginary part of the perturbed lowest eigenvalue. 
Dotted line: perturbation data. Solid line: data from [g. 



a numerical integration of a differential equation that 
presents some accuracy problem for very low values of 
the coupling constant and works better at higher values: 
the imaginary part is now deduced from the half- width 
of the Breit-Wigner lines that fit the continuous spec- 
trum. This makes the two methods complementary and 
puts rather narrow bounds to the range of the f2 val- 
ues in which the comparison makes sense: a reasonable 
interval for the present data could be assumed to be 
0.1 < fi < 0.25. Within these bounds Fig.l shows that 
the agreement is complete, proving the effectiveness of 
both the DBS and the method of [y] for solving numer- 
ical problems. 



APPENDIX A. 



r) 

|Re(a 2 



(Im(o-)y 2 + ilm(a)Im(a 2 ) y 5 + J^y 8 ) 



.14 



16 



(p 2 + y 2 + Re(a)y±^y 4 ) 2 



+ Im(cr)p + Im(cr 2 )(py 3 + y s p) 



For < a< l—\Re(a 2 )\/\a\ 2 and for a certain R> we 
then have 



(To + Ti)* (To + Ti) > a [(p 2 + y 2 + Re(a)y) 2 

+ (Im(a)y 2 + ~Im(a)Im(a 2 )y 5 + J^y 8 ) 
4 16 

' |Re(fj2)l [lm(a)y 2 + ^lm(a)lm(a 2 )y 5 



R 

Im(cr 



(p 2 + y 2 + Rc(o-)y) 2 + ^Ly* 
■ (P± 1) 2 -^(P 2 + 1) 



2 ^ ' 2 
\lm(a 2 )\(p±y 3 ) 2 -\lm(cr 2 )\(p 2 + y e ) 



> a (T* T + T* Tj-b+R ( P 2 + y 2 + Re(a)y) 



|Im( CT 2 ) b 2 -^ 2 



|Re(a 2 )| /l 



r |4 



(-Im(a)Im( ( 7 2 )y 5 + Im((j)y 2 ) 

b 
2 



R^-y*-\l^)\y« 



(A.l) 



In this Appendix we give the proof of the quadratic 
inequality pil.ip with T and T x given in pn.6[) . We 
have: 

(To + Ti)* (T + Tl) = (p 2 + y 2 + Re(a)y) 2 

+ (* Im(cr)y + -j2/ 4 ) Im ( cr )y + — V A ) 

2 

+ (P 2 + y 2 + Re(a)y) («Im(a)y + ^y 4 ) 

- 2 

+ (-* Jm(a)y + ^y 4 ) (p 2 + y 2 + Re(a)y) 



In the last inequality we have neglected two positive 
terms and we have chosen a number b > in such a way 
to make positive the last two square brackets. In spite 
of the explicit calculation, one could also have argued 
that the linear term ay is small in the sense of Kato 
with respect to the leading term of the perturbation 
(<r 2 /4)y : consequently it gives a negligible contribu- 
tion to the quadratic estimate which therefore reduces 
to that of the anharmonic oscillator. We have chosen 
this shortcut to deal with the quadratic estimates of the 
three-dimensional case. 
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